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Abstract 

Using the concept of open systems where the classical geometry is treated as the 
system and the quantum matter field as the environment, we derive a fluctuation- 
dissipation theorem for semiclassical cosmology. This theorem which exists under very 
general conditions for dissipations in the dynamics of the system, and the noise and 
fluctuations in the environment, can be traced to the formal mathematical relation 
between the dissipation and noise kernels of the influence functional depicting the 
open system, and is ultimately a consequence of the unitarity of the closed system. In 
particular, for semiclassical gravity, it embodies the backreaction effect of matter fields 
on the dynamics of spacetime. The backreaction equation derivable from the influence 
action is in the form of a Einstein-Langevin equation. It contains a dissipative term 
in the equation of motion for the dynamics of spacetime and a noise term related to 
the fluctuations of particle creation in the matter field. Using the well-studied model 
of a quantum scalar field in a Bianchi Type-I universe we illustrate how this Langevin 
equation and the noise term are derived and show how the creation of particles and 
the dissipation of anisotropy during the expansion of the universe can be understood 
as a manifestation of this fluctuation-dissipation relation. 
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1 Introduction 



An important relation between the dissipation in the dynamics of a system and the fluctu- 
ations in a heat bath with which the system interacts is the fluctuation-dissipation relation 
(FDR) |J. A first example of its manifestation is the Nyquist noise in an electric circuit ||. 
This relation is of practical interest in the design of noisy systems ||. It is also of theoreti- 
cal interest in statistical physics because it is a categorical relation which exists between the 
stochastic behavior of many microscopic particles and the deterministic behavior of a macro- 
scopic system. It is therefore also useful for the description of the interaction of a system 
with fields, such as effects related to radiation reaction and vacuum fluctuations between 
atoms and fields in quantum optics [|J. The form of the FDR is usually given under near- 
equilibrium conditions via linear response theory (LRT) ||. We will see in this paper that 
this relation has a much wider scope and a broader implication than has been understood 
before. In particular we want to address problems involving gravity and quantum fields in 
black holes and the early universe; and we are interested in seeing this relation validated and 
implemented under non-equilibrium conditions for quantum fields in curved spacetimes 0. 
The problem we choose for illustration is the backreaction of particles created J7|, |8|, [10], [11 
in a cosmological spacetime (without event horizon) |L^, [13], [14], [15], [L6|, [T7[ [L8|, [LJJ]. The con 



ceptual framework we adopt is that of a quantum open system [2C], the formal scheme is 
that of the closed-time-path and the influence functional [22 , formalisms, and the 



paradigm we use for comparison is that of quantum Brownian motion [|24| , |25 
as the one who had the great insight of proposing a flucti 
for the depiction of quantum processes in black holes 



Sciama [27| was the one who had the great insight of proposing a fluctuation-dissipation 
relation |28, 

if and de Sitter universe 



accelerated observers 31 



OJ , uniformly- 
Using Einstein's analysis of the 
Brownian motion as a guide he showed that the Hawking and Unruh radiations can be seen 
as excitations of vacuum fluctuations and the detector response as following a dissipation- 
fluctuation relation. A crucial element for this interpretation to be possible is the existence 
of an Euclidean section in the Schwarzschild, Rindler and de Sitter metrics, imparting a 
periodicity in the Green's function of the matter field, thus turning it into a thermal prop- 
agator 0, [35[ (in the imaginary-time Matsubara sense), and forging the equivalence of the 
system with finite-temperature results P5| . The same condition applies to Hawking radi- 
ance in de Sitter universe [[HJ, which by virtue of its possession of an event horizon, also 
admits a FDR interpretation f37| . The derivation of the FDR in this class of spacetime was 
based on a linear response theory, which hinges on the thermal equilibrium condition set 
up by the created particles. For spacetimes without an event horizon, or for systems under 
non-equilibrium conditions, one would not ordinarily think that a FDR could exist |5SL ^UJ . 
The generalization of this relation to non-equilibrium conditions is a much more difficult 
problem. 

When Sciama first proposed this way of thinking, one of us was involved in the backreac- 
tion studies of quantum processes in cosmological spacetimes [[HJ. The dissipative effect of 
particle creation on the dynamics of spacetime seems to point to the existence of a fluctuation- 
dissipation relation, except that one factor (dissipation) is not clear, and the other factor 
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(noise) is missing. Two important steps had to be taken before this picture began to make 
better sense. In the calculation of Hartle and Hu |17| on anisotropy damping, the Schwinger- 
DeWitt effective action ('in-out') formalism |H| gives rise to an effective geometry which is 
complex, making it difficult to interpret what dissipation really is. The adoption of the 
Schwinger-Keldysh (closed-time-path, CTP, or 'in- in') [^TJ formalism by Calzetta and Hu 
iQfl yields a real and causal equation of motion for the effective geometry, from which one 



can relate the source of dissipation in spacetime dynamics to the energy density of particles 
created and explicitly identify the viscosity function associated with anisotropy damping 



41 



The adoption of the CTP formalism was an encouraging step in the right direction, but 
one needs to understand the statistical mechanical meaning of these quantum processes bet- 
ter in order to appraise the validity of adopting the well-established concepts and results in 
statistical mechanics for their depiction or explication. For the particular task of showing a 
fluctuation-dissipation relation at work for quantum fields in a general cosmological space- 
time not required to possess an event horizon, there is also the noise or fluctuation term 
missing. These inquiries were summarized in a report written by one of us [4~2 in which 
some tentative replies were given in the form of three conjectures: 

1) That colored noise associated with quantum field fluctuations is generally expected in 
gravitation and cosmology; 

2) That the backreaction of particle creation in a dynamical spacetime can be viewed as the 
manifestation of a generalized fluctuation-dissipation relation; and 

3) That all effective field theories, including semiclassical gravity or even quantum gravity 
(to the extent that it could be viewed as an effective field theory), are intrinsically dissipative 
in nature. 

There, it was also suggested that one can use the Caldeira-Leggett model |Z3 to study the 
theoretical meaning of dissipation and probe into the relation of noise and dissipation. The 
next stage of work in this quest concentrated on the properties of quantum open systems 
PDfl and extending the theory to quantum fields and to curved spacetimes. 

Using the influence functional formalism of Feynman and Vernon p2 |, one can identify 
noise in an environment from the imaginary part of the influence action. The characteristics 
of noise depends on the spectral density of the environment, the coupling of the system 
with the environment, and other factors. Using a model of the Brownian particle coupled 
nonlinearly with a bath of harmonic oscillators, Hu, Paz and Zhang [25| deduced the noise 
autocorrelation functions and a generalized fluctuation-dissipation relation for systems driven 
by intrinsic colored and multiplicative noises. This forms the basis for the second stage of 
investigation. To generalize these results to quantum fields, Hu and Matacz [43[ 44| recently 
analyzed the problem of QBM in a parametric oscillator bath. A parametric oscillator 
bath enables one to study particle creation in quantum fields, where the Brownian particle 
can play the role of an Unruh detector, or, in a cosmological backreaction problem, the 
scale factor of the universe. One can study the detector's response or the effect on the 
universe due to the fluctuations of the quantum field. They found that the characteristics 
of quantum noise vary with the nature of the field, the type of coupling between the field 
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and the background spacetime, and the time-dependence of the scale factor of the universe. 
They showed how a uniformly accelerating detector in Minkowski space, a static detector 
outside a black hole and a comoving observer in a de Sitter universe all observe a thermal 
spectrum. By writing the influence functional in terms of the Bogolubov coefficients which 
determine the amount of particles produced, they also identified the origin of noise in this 
system to particle creation |4"5| , |46| . The influence functional method not only reproduces 
the known results, but also enables one to look into the hitherto unknown domain of noise, 
fluctuations, and decoherence. 

A program for studying the backreaction of particle creation in semiclassical cosmology 
in the open system conceptual framework using influence functional methods was recently 



outlined in |46|, |J, [47[]. The backreaction of these quantum field processes manifests as 
dissipation effect, which is described by the dissipation kernel in the influence action. Using 
a model where the quantum Brownian particle and the oscillator bath are coupled paramet- 
rically (the field parameters change in time through the time- dependence of the scale factor 
of the universe, which is governed by the semiclassical Einstein equation) Hu and Matacz 
f4"S| derived an expression for the influence functional in terms of the Bogolubov coefficients 
as a function of the scale factor. From the variation of the influence action they obtained an 
equation of motion describing the dynamics of spacetime in the form of an Einstien-Langevin 
equation. 

After these recent works, it is clear that the influence functional method is the appropriate 
framework for studying the nature and origin of noise in quantum fields and to explore the 
statistical mechanical meaning of quantum processes like particle creation and backreaction 
in the early universe and black holes. Two additional aspects, however, need be considered 
to complete the story. First, how is it related to the CTP formalism, which gave us, to begin 
with, the correct dissipation side of the story? This problem was taken up in a recent paper 
of Calzetta and Hu [[45], who showed how noise and fluctuations in semiclassical gravity 
can also be obtained with the original CTP formalism. They also showed that the CTP 
and the IF formalisms are indeed intimately related. They derived an expression for the 
CTP effective action in terms of the Bogolubov coefficients and showed how noise is related 
to the fluctuations in particle number. From there, they show how an Einstein-Langevin 
equation naturally arises as the equation of motion for the effective geometry, from which 
a new, extended theory of semiclassical gravity is obtained. The work of Calzetta, Matacz 
and the present authors shows clearly that the old framework of semiclassical gravity is 
only a mean field theory. This theory based on the Einstein equation with a source driven 
by the expectation value of the energy momentum tensor should be replaced by one based 
on an Einstein-Langevin equation which describes also the fluctuations of matter fields and 
dissipative dynamics of spacetime. 

Notice that in moving from the first stage of this investigation based on the CTP for- 
malism to the second stage based on the IF formalism, one has to elevate the treatment of 
classical spacetimes as external fields to reduced density matrices. In making these transi- 
tions and back, several issues need be addressed. The central issue is the quantum to classical 



transition for the spacetime sector 4£ . The important question behind the transition from 
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quantum gravity to semi classical gravity is decoherence. This is a subject of much recent 
interest. We refer the reader to recent work for the exposition of different viewpoints and 
approaches K9], [50L RTF] . Here, for the backreaction problem, we shall adopt the results of 



Paz and Sinha [ 52]j , which is based on a reduced density matrix formalism adapted to quan- 
tum cosmology. There, the model of a Bianchi-I universe coupled to a scalar field was used 
to derive conditions for transition from quantum cosmology to the semiclassical limit via 
decoherence, and the relationship between decoherence and backreaction was investigated. 

After these previous investigations paved the way for the use of open-system concepts 
applied to the backreaction problem of quantum fields in curved spacetimes, we are finally 
in a position to look at the full picture and explore the existence of a fluctuation-dissipation 
relation for semiclassical gravity in general. We shall use the model of particle creation 
in Bianchi Type I universe to explore this relation. In Sec. 2, we give a summary of the 
results for the quantum Brownian model, assuming a general nonlinear coupling between 
the system and the environment, giving rise to colored and multiplicative noise. Readers 
familiar with the QBM problem can skip over this section. In Sec. 3 we begin with the 
density matrix of the universe and show how coarse-graining the matter field viewed as an 
environment produces the reduced density matrix, and how the influence functional defined in 
the evolutionary operator for the reduced density matrix contains the relevant information 
we need- the dissipation and noise kernels. In Sec. 4 we analyze the phase and the real 
components of the influence functional in detail, sorting out the divergent and renormalized 
terms in the phase. We show that the renormalized phase part provides the dissipative term 
in the equation of motion, and the real component contributes to decoherence and noise. We 
show how a colored noise of the quantum field can be identified, and with it the existence of 
a fluctuation-dissipation relation between these kernels. In Sec. 5, we discuss the physical 
meaning of this relation. We first show that noise measures the difference in the amounts 
of particle creation along two histories. Since this is also the condition for decoherence to 
occur, we see that a relation also exists between decoherence and particle creation. With 
this noise term, we then derive the Einstein-Langevin equation for the anisotropy tensor. We 
show that it is identical in form to that derived via the CTP formalism before |4(|, but with 
a new stochastic source term from the noise, as anticipated in Finally, we show how 
the dissipation in the anisotropy of spacetime can be related to the particles created. Thus 
noise and dissipation which are connected by a formal relation, are both related to particle 
creation, and the backreaction of particle creation is an embodiment of the FDR. In Sec. 6 
we discuss the physical interpretation of the FDR in a more general context. We show how 
the changing rate of particle creation and the strength of backreaction effect can be gauged 
consistently by the fluctuation-dissipation relation valid for time- dependent conditions. We 
also describe related problems for future investigations. 



2 Influence Functional for Quantum Open System 
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2.1 Quantum Brownian Motion Paradigm 

Let us first review a model problem of quantum Brownian motion (QBM) where the role of 
noise and dissipation are well understood. Subsequently we will draw analogies from this 
problem to analyze the quantum cosmology problem of our interest. 

Consider a Brownian particle interacting with a set of harmonic oscillators. The classical 
action of the Brownian particle is given by 

S[x] = J* ds{^Mx 2 -V{x)}. (2.1) 

The action for the environment is given by 

Se[{q n }} = J q dsY,{^rn n q 2 n - ^m n ulql}. (2.2) 



We will assume that the action for the system-environment interaction has the following 
form 

t 

Sint[x,{q n }} = / ds^v n (x)q k n (2.3) 
o 

where v n (x) = —Xc n f(x) and A is a dimensionless coupling constant. If one is interested 
only in the averaged effect of the environment on the system the appropriate object to study 
is the reduced density matrix of the system p r , which is related to the full density matrix p 
as follows 



-oo +oo 



p r (x,x') = J dq J dq'p{x,q;x',q')5{q-q'). (2.4) 

— oo — oo 

It is propagated in time by the evolution operator J r 

+oo +oo 

p r (x,X,t)= J dXi J dx[ J r (x, X , t | Xi, x[, 0) /0 r (Xj, x[, ). (2.5) 

— oo — oo 

If we assume that at a given time t — the system and the environment are uncorrelated 

p(0) = p s (0)xp e (0), (2.6) 
then J r does not depend on the initial state of the system and can be written as 

x f x 'f 

J r {xf,x'j,t | Xi,x\,) = J Dx J Dx' exp— js*^] — S^a/]} J-[x, x'\ 



x' 



j Dx j Dx' exp^S eff [x,x'} (2.7) 
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where the subscripts i, f denote initial and final variables, and S e ff[x,x'] is the effective 
action for the open quantum system. The influence functional J-[x, x'} is defined as 



+00 +00 +00 if if 
T\x,x'\ = J dqj J dqi J dq[ J Dq J Dq' 



x exp ^{S e [q\ + S int [x, q\ - S e [q] - S int [x' ', q']}p e (qi, ?*, 0) 



exp — Sif[x, x'\ 



{21 



where Sjp[x,x'} is the influence action. Thus S e ff[x, x'] = S[x] — S[x'] + Sjp[x,x']. 

From its definition it is obvious that if the interaction term is zero, the influence functional 
is equal to unity and the influence action is zero. In general, the influence functional is a 
highly non-local object. Not only does it depend on the time history, but -and this is 
the more important property- it also irreducibly mixes the two sets of histories in the path 
integral of (2.7). Note that the histories x and x' could be interpreted as moving forward and 
backward in time respectively. Viewed in this way, one can see the similarity of the influence 
functional and the generating functional in the closed-time-path, or Schwinger-Keldysh pi[ 
integral formalism. 

We will assume that initially the bath is in thermal equilibrium at a temperature T = 
{k B (3)~. The T = case corresponds to the bath oscillators being in their respective ground 
states. It can be shown |26| that the influence action for the model given by the interaction 
in (2.3) to second order in A is given by 



S IF [x,x'} = {J ds [-AV(x) }- J ds [-AV(x')} } 



fc HI 



dso A^ 



f(x( Sl ))-f(x'( Sl )) 



(*), 



si - s 2 ) 



f(x(s 2 )) + f(x'(s 2 )) 



+i / dsi / ds 2 A 



(*), 



si - s 2 



f(x(s 2 ))-f(x>(s 2 )) (2.9) 



f( x ( Sl ))-f(x'( Sl ))\u 


where AV(x) is a renormalization of the potential that arises from the contribution of the 
bath. It appears only for even k couplings. For the case k = 1 the above result is exact. 



This is a generalization of the result obtained in p2 where it was shown that the non-local 



kernel pS '{s\ — s 2 ) is associated with dissipation or the generalized viscosity function that 
appears in the corresponding Langevin equation and u^ k \s\ — s 2 ) is associated with the 
time correlation function of the stochastic noise term. The dissipation part has been studied 
in detail by Calzetta, Hu, Paz, Sinha and others 



45 1 in cosmo logical 
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backreaction problems. We shall elaborate somewhat on the nature of the noise part of the 
problem and then analyze their connection. In general v is nonlocal, which gives rise to 
colored noises. Only at high temperatures would the noise kernel become a delta function, 
which corresponds to a white noise source. Let us first see the meaning of the noise kernel. 
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2.2 Noise 

The noise part of the influence functional is given by 



exp{-- / ds, I ds 2 \f(x( Sl )) - f(x'( Sl ))y k \ Sl - s 2 )\f(x(s 2 )) - f(x'(s 2 ))} (2.10) 







where is redefined by absorbing the A 2 . This term can be rewritten using the following 
functional Gaussian identity [^] which states that the above expression is equal to 



. t 

J V^\t)V[^]eM^Jds^ k \s)[f(x(s))-f(x f (s))} 



(2.11) 



where 



1 * * 1 

p[£W] =P«exp{-~| d Sl J ds 2 -^ k \s 1 )[^ k \s 1 -s 2 )]~ 1 ^ k \s 2 )} (2.12) 







is the functional distribution of ^ k \s) and is a normalization factor given by 



t t 



1 

[p(*>]-i = | V^ k \s)exp{--J ds x f ds 2 ^ k \ Sl )[v^{ Sl - s 2 )r l ^ k \s 2 )}. 



(2.13) 





The influence functional can then be rewritten as 



2?eW(a)P[e ( *Wp/F[a:.ar'.e C * ) ] 



exp^^z',^ 



(2.14) 



where 



t t 
S IF [x,x',^ k) ] = J ds {-AV(x) } - I ds {-AV{x') } 



t Si 

- jds x J ds 2 [f(x( Sl )) - f(x'( Sl ))}^ k \ Sl - s 2 ) [f(x(s 2 )) + f(x'(s 2 )) 



t t 



-j dse k \s)f(x(s)) + j ds^ k \s)f(x'( S )) 



so that the reduced density matrix can be rewritten as 

Pr (x,x') = J V^ k \s)V[^] Pr (x,x',[^}). 



(2.15) 



(2.16) 



The full effective action can be written as 

t 



S eff [x,x',Z] = {S[x} + / ds Z(8)f{x{s))} - {S[x'\ + / ds i{s)f{x'{s))} 





si 

ds x Jds 2 [f{x{s x )) - f( x '(s 1 ))y k \s 1 - s 2 )[f{x(s 2 )) + f(x'(s 2 $]l7) 



From equation (2.15) we can view £^(s) as a nonlinear external stochastic force and 
the reduced density matrix is calculated by taking a stochastic average over the distribution 
.P[£^] of this source. 

From (2.12), we can see that the distribution functional is Gaussian. The Gaussian noise 
is therefore completely characterized by 



(2.18) 



We see that the non-local kernel i/W(si — s 2 ) is just the two-point time correlation function 
of the external stochastic source C ( s ) multiplied by %. 

In this framework, the expectation value of any quantum mechanical variable Q(x) is 
given by |51 



(Q(x)) = Jv^ k \s)V[^] J dx Pr (x,x,[e ] ])Q(x) 

— oo 

'M^uantum)^ (2-19) 

This summarizes the interpretation of u^(si — s 2 ) as a noise or fluctuation kernel. 



2.3 Langevin Equation 

We now derive the semiclassical equation of motion generated by the influence action (2.9). 
This will allow us to see why the kernel ^ k '(sx — s 2 ) should be associated with dissipation. 
Define a "center-of-mass" coordinate x and a "relative" coordinate A as follows 



x[s) 

Ms) 



~[x(s)+x'(s)} 



- x'(s) — x(s). 

The semiclassical equation of motion for x is derived by demanding (cf. M 

5 r 



5A 



S[x] — S[x'\ + Sjf[x,x') 



0. 



A=0 



(2.20) 



(2.21) 
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Using the sum and difference coordinates (2.20) and the influence action (2.9) we find that 
(2.21) leads to 

dL r ddL r df(x) }, a (k) ( . ^ df(x(s)) 



where 4- r y^{t — s) = fj,W(t — s). We see that this is in the form of a classical Langevin 

equation with a nonlinear stochastic force F^k)(s) = —^ k \s) 9 ^ . This corresponds to a 
multiplicative noise unless f(x) = x in which case it is additive. L r denotes a renormalized 
system Lagrangian. This is obtained by absorbing a surface term and the potential renor- 
malization in the influence action into the system action. The nonlocal kernel j^(t — s) is 
responsible for non-local dissipation. In special cases like a high temperature ohmic environ- 
ment, this kernel becomes a delta function and hence the dissipation is local. 



2.4 Fluctuation-Dissipation Relation 

Recall that the label k is the order of the bath variable to which the system variable is 
coupled. 7 < - fc - ) (s) can be written as a sum of various contributions 

7 (fc) (*)=£7z (fc) (s) (2-23) 
i 

where the sum is over even (odd) values of / when k is even (odd). To derive the explicit 
forms of each dissipation kernel, it is useful to define first the spectral density functions 

= £ % - u n ) k *h™ p^- k - (2.24) 

It contains the information about the environmental mode density and coupling strength 
as a function of frequency. Different environments are classified according to the functional 
form of the spectral density I(oS). 

In terms of these functions, the dissipation kernels can be written as 

(fc) = r au i j(fc) ^ M {k)^ z) cQslus (2 _ 25) 



where Mi k \z) are temperature dependent factors derived in [26| and z = coth^(3hu. Anal- 
ogously, the noise kernels u^ k \s) can also be written as a sum of various contributions 

^(sHE^OO (2-26) 
I 

where the sum runs again over even (odd) values of I for k even (odd). The kernels v x f (s) 
can be written as 



.(*) 



+oo 

h [ — lW(u) N} k \z) cos Iujs (2.27) 

J IX 
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where (z) is another set of temperature- dependent factors given by 

To understand the physical meaning of the noise kernels of different orders, we can think 
of them as being associated with I independent stochastic sources that are coupled to the 
Brownian particle through interaction terms of the form (2.15) 



fdsJ2 /(*)■ (2-28) 

o 1 



This type of coupling generates a stochastic force in the associated Langevin equation 

F^(s) = -tf\s)^ (2.29) 

which corresponds to multiplicative noise. The stochastic sources £j have a probability 
distribution given by fl2.12|) which generates the correlation functions ( |2.18| ) for each k and 
I. 

To every stochastic source we can associate a dissipative term that is present in the real 
part of the influence action. The dissipative and the noise kernels are related by generalized 
fluctuation-dissipation relations of the following form 

+ 0O 

„«>•(*)= jdsK^if - S)7 ?' W (2.30) 

— OO 

where the kernel k\ (s) is 

+oo 

K\ k \s) = [ — L?\z) I u cos Iujs (2.31) 

J 7T 



and the temperature-dependent factor L\ k \z) = N' t (z) / 'Mi k \z) . 

A fluctuation dissipation relation of the form ( |2.30| ) exists for the linear case where the 
temperature dependent factor appearing in (|2.31|) is simply = z. The fluctuation- 
dissipation kernels have rather complicated forms except in some special cases. In the 
high temperature limit, which is characterized by the condition k^T ^> hA, where A is the 
cutoff frequency of the environment, z = coth f3%uj / '2 — > 2/(3Hu we obtain 

L?\z) - • (2.32) 
In the limit A — > +oo, we get the general result 

K\ k \s) = 2 -^fs(s) (2.33) 
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which tells us that at high temperature there is only one form of fluctuation-dissipation 
relation, the Green-Kubo relation Jl| 

4%) = ^>>M. (2.34) 

In the zero temperature limit, characterized by z — > 1, we have 

Lf\z)^l. (2.35) 

The fluctuation-dissipation kernel becomes ^-independent and hence identical to the one for 
the linearly- coupled case 

+ °° d 

K( s ) = / — ujcoslus. (2.36) 

J IX 


It is interesting to note that the fluctuation-dissipation relations for the linear and the 
nonlinear dissipation models are exactly identical both in the high temperature and in the 
zero temperature limits. In other words, they are not very sensitive to the different system- 
bath couplings at both high and zero temperature limits. The fluctuation-dissipation relation 
reflects a categorical relation (backreaction) between the stochastic stimulation (fluctuation- 
noise) of the environment and the averaged response of a system (dissipation) which has a 
much deeper and universal meaning than that manifested in specific cases or under special 
conditions. 

Our aim in the next section would be to consider a model consisting of quantum fields 
coupled to a cosmological background metric and cast it into the system-environment form 
as discussed here. Consequently we shall see that one can construct an influence functional 
of a form very similar to (|2.8|) and hence derive a fluctuation-dissipation relation of the form 



3 Influence Functional for Quantum Cosmology 
3.1 Reduced Density Matrix of the Universe 

The model we will analyze here is the same as that used in [f}2] from which we will quote 



results relevant to our study. Our "system" will consist of a minisuperspace model with D 
degrees of freedom denoted by coordinates r m (with m = 1, . . . ,D). The minisuperspace 
modes will be coupled to "environment" degrees of freedom that we schematically represent 
by <3> (they will be later associated with the modes of a scalar field) . The quantum mechanical 
description of this Universe will be given by the wave function of the Universe ^ = \l/(r m , $) 
which, as a consequence of the existence of a classical Hamiltonian constraint, satisfies the 
Wheeler- DeWitt equation: 

H^I = (H r + H$)ty = (3.1) 
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In the class of models we consider, the Hamiltonian corresponding to the minisuperspace 
variables can be written as 

H r = ^G mm ' PmPm , + MV(r m ) (3.2) 

The matrix G mm ' determines the metric in the minisuperspace (the supermetric) and the 
quantity M is proportional to the square of the Planck mass. In the following we will set 
h = 1 throughout. In the above Wheeler- DeWitt equation we assume that the momenta are 
replaced by operators according to a covariant factor ordering prescription. The Hamiltonian 
constraint represents an important distinction from the quantum Brownian motion case 
discussed previously, because it implies that there is no preferred notion of time in this case 
and the wavefunction satsfies (|3.1|) rather than the Schrodinger equation. The Hamiltonian 
associated with the environment degrees of freedom is some function -£/$($, 7r$, r m ,p m ) that 
we will specify later. 

We will be interested in making predictions concerning only the behavior of the minisu- 
perspace variables r m which we consider the "relevant" part of the universe. To achieve such 
a coarse- grained description we will work with the reduced density matrix of the system 
which is defined as: 

Pred(r', r) = J d$**(r, $)*(/, $) (3.3) 

For some region of the minisuperspace, (|3.1|) admits solutions that are oscillatory functions 
of r m of the following WKB form: 

ty(r, $) = e iMS{r) C{r)^j{r, $) (3.4) 

In this regime, the system variables r and the environment variables <3? behave as heavy 
and light modes respectively (the Planck mass plays the role of a large mass parameter) in 
analogy with the Born-Oppenheimer approximation. This also provides some justification of 
the system-environment split akin to the Brownian motion case . Thus if one assumes that 
all the functions S, C, ip can be expanded in powers of M _1 and substitutes these expansions 
into (|3.1|), one gets, to leading order (i.e., M°): 

^ m " , '^^ + ^M = (3.5) 

which is essentially the minisuperspace version of the Hamilton- Jacobi equation. To the next 
order in M one obtains, 

iGmm '^^ = = (3.6) 

This last equation is obtained provided we choose the prefactor Co identical to the H$ = 
case. Thus, if we define the WKB time t as 

j _ Qmm' dffp 9 . . 

dt dr m ' dr m 1 ' ' 
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the equation ( |3.6| ) reduces to the familiar Schrodinger equation that reads: 

l lt = H ^ (3 ' 8) 

From now on we will drop all the O-subindices which should be considered as implicit in all 
the equations where S, C and if) appear. The Hamilton- Jacobi equation ( |3.5| ) will have a 
D — 1 parameter family of solutions and for each one of these solutions we can build a wave 
function like (3.4|). In general one can assume that the wave function of the Universe is a 



superposition of these terms, each of which will be called a WKB branch: 

*(r,<t>) ^ ei%)(r)C w(^(»)M) ( 3 - 9 ) 

n 

Here the subindex (n) labels the WKB branch characterized by a set of parameters (n) that 
uniquely defines the particular solution to the Hamilton- Jacobi equation. However, in the 
rest of our analysis , we will consider the wavefunction to be represented by a single term of 
the above sum, i.e, by a particular WKB branch. We will drop the subscript n from now on 
with this understanding. 

The reduced density matrix associated with the wave function (|3.4j) is: 



Pred{r', r) = e iM[s{r) - s{r ' )] C(r)C(r')l(r', r) (3.10) 

where 



l(r',r) = J ^*(r',$)V<r,$)d$ (3.11) 
The influence of the environment on the system is summarized by the above function X and 



it will be the basic object of our interest. It has been shown in references p2| , [55j that this is 
the object that is exactly analogous to the Feynman- Vernon influence functional T[x^ x') in 
the case where the environment is in a pure state. We will therefore call X(r', r) the influence 
functional and analyze the fluctuation and dissipation phenomena in analogy to the QBM 
problem. To facilitate making these connections, we write the influence functional in the 
form 

J(r, r') = exp{ir(r, r')} (3.12) 
where the influence action can be written as 

r(r,r) = 6(r,r') + it(r,r') (3.13) 

The phase and the real exponent T which constitute the influence functional will be the 
basic objects of our interest (note that T is positive since the overlap is bounded by unity). 

3.2 Bianchi-I Minisuperspace with a Conformal Scalar Field 

We now specialize our model to a minisuperspace of Bianchi I universe coupled to a massless 
conformal scalar field. The line element is given by []5f| 

ds 2 = a 2 dr] 2 - c^ef.dx'dx^ (3.14) 
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as: 



where r] is the conformal time . The traceless 3x3 matrix (3 measures the anisotropy, its 
time rate of change gives the shear. For Type-I universe, it can always be parametrized by 
the principal eigenvalues 

/3 = diag(/5 1 ,/3 2 ,/3 3 ) (3.15) 

or, equivalently by (3± defined by 

A = /? + + V3/3_, fa = (3 + - y/3/3-, /3 3 = -2/5+ (3.16) 
Rewriting the scale factor as a = e a , the Einstein Hilbert action can be written as 

S g = 6M J dr]{e 2a {-a 2 + + $ 2 _) (3.17) 

where M = Mp, , and a dot denotes taking a derivative with respect to the conformal time 
f). We normalize the spatial volume to 1 assuming T 3 spatial topology. 
The action for the scalar field is given by 

Sf = I J d 4 x (g^d^d^ - ii?$ 2 ) (3.18) 
which, after integrating by parts and defining the conformal field X = a$, can be written 

S f = ^J d 4 x {X 2 + XV (3) X - 01 + (3 2 _)X 2 } (3.19) 

where the spatial Laplacian is given by = e^didj. 

As usual, we expand the field X in an orthonormal basis of eigenfunctions of V^. As 
the spatial sections are flat, the eigenfunctions are simple trigonometric functions and the 
momenta are quantized due to the periodic boundary conditions associated with the T 3 . We 
will denote the basis as {Qko-(af), k = (k x , k y , k z ), kj = 2irnj, a = ±}. The index a labels the 
functions according to their parity. The expansion of the field X reads: 

X(x, V ) = J2QUZ)xUv) (3-20) 

kcr 

The variables of the minisuperspace constituting our open system are r m = (a, /3 + ,/3_) 
or (a,Pij) and the 'environment' variables are the collection of field amplitudes {xka, kj = 
27rrij, a = ±}. Using our previous expressions it is easy to show that the Hamiltonian can 
be written in the form of (|3.2| ), where the gravitational part has the supermetric 

<w = ldiag(-l,+l,+l). (3.21) 

or 

On the other hand the matter Hamiltonian can be written as 

H x = £^ = E^L + ^Xk CT 2 ) 

kcr kcr 

^ = +Tu^ {p ^ +p i- ) (3 ' 22) 
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We will assume the wave function of the universe can be written as ( |3.4| ), where the 
function S obeys the Hamilton-Jacobi equation (|3.5|) which in this case is given by: 



e -2a 



(-(d a S) 2 + (dp + S) 2 + (dp_Sf) = (3.23) 



2M 

This equation can be separated and solved as 

S(a, p±) = S % {a) + b + (3 + + b_{3_ (3.24) 

with 

dJzia) = ±\b\ (3.25) 

where we use b to denote the two dimensional constant vector (6 + ,6_). 

As we can see, a particular solution to the Hamilton-Jacobi equation is parametrized 
by two integration constants (b + and &_) and by the sign that defines S(a) in equation 
( |3.24j ). Therefore, the label (n) that characterizes a solution of ( |3.23j ) stands for the set of 
constants {b, ±}. Every function S( n ) generates a 2-parameter family of trajectories in the 
three dimensional minisuperspace (these are the curves orthogonal to the S7„) =constant 
hypersurfaces) . These trajectories are exact solutions to the Einstein's equations, and if we 
restrict our considerations to the plane (a, (3 + ), the trajectories are defined by the equation 

^± - — b -+- (3 26) 

da ~ d a S {6 - Ib) 

The minisuperspace trajectories can be found by integrating the above equation and are 
straight lines (with slope given by ±6 + /|6|) corresponding to the well known Kasner's solu- 
tions. In that case, for the "expanding" (i.e. a > 0) branch, we have (3 + = + /3 +0 , where 
f3 +0 is an integration constant. 



3.3 Influence Action 



We have to compute the influence functional ( |3.13| ) according to the strategy described in 
the beginning of this section and for that we have to solve the Schrodinger equation ( |3.6| ). 
It is possible to make the following ansatz for the matter wave function 

V(r, X) = V(r, {Xk}) = II Mr, Xk) (3-27) 

k 

Thus, the influence functional is expressed as an infinite product while the phase and the 
real exponent T can be written as a sum of contributions from each mode. 

Each component of the wave function satisfies the following Schrodinger equation: 

i^r = H * ^- ( 3 - 28 ) 

or] 
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with a Hamiltonian given by 



= -~ + \{e^ M + fc)xl (3.29) 

where as before, we have used a dot to denote the derivative with respect to the conformal 
time, which also happens to coincide wth the WKB time as can be seen from applying the 
definition (3.7) to the model of sect. 3.2. 

Let us now describe how we compute the influence functional. We will make a Gaussian 
ansatz for the wave function ip^ that corresponds to assuming that the state for the scalar 
perturbations is a particular vacuum . Thus, we write each component of the wave function 
as (for simplicity we will omit the index k): 

^rj) = (-)i e -U^ dt e ^ 2w (3.30) 

where w = ii/u = w r + iwi, and w r , Wi are the real and imaginary parts of w. The equation 
satisfied by the function u is easily derived from the Schrodinger equation and can be written 
as : 

u + ttlu = (3.31) 

The computation of the overlap factor involves solving the above equations. In our model 
this can be done using a perturbative scheme if we assume that the anisotropy coordinates 
are small. In that case, we can can write (up to second order in the anisotropy): 



where 



Qi = ui - Ai - A 2 (3.32) 



|k 2 | 1/2 , Ai = -2/%jW and 



A 2 = - 2Plk t h? -{Pl + 0l) (3.33) 



Then, the equation for u can be solved by a standard iteration procedure |52|, |3T|. Assuming 
that the anisotropy is "switched off" at early and late times, and taking the initial state as the 
conformal vacuum, the expressions for T and of the exponent of the influence functional 
defined in fl3.13Q are given respectively by 



T(r,r') = ulcos (2u; k - r/ 2 )) 

1 r' j Ai(*7i)Ai(tfc) 



lb J J out 



~ r f dm d V2 Xlivi) f ir]2) cos(2u k ( v ' -v + m- m)) (3.34) 

lo J J out 
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and 



1 1 rv In' 

®(r,r') = ~u> k (r] - r\ ) + - — / dr\x\ 2 {T]x) - - — / dr}\\ 2 (T}x) + 



2 4cu k J Auj v 

1 rr, rvi A X (?7i) \i{r) 2 



+ - / / dr]i dr] 2 5 szn(2a; k (^i - r] 2 )) ~ 

8 J J UJy. 

- - / dr]i dr] 2 5 sm(2a> k (77i - 772)) + 

o J J to> k 

+ i /" r *n * h -(2^(V - , + * - *)) (3.35) 



up to second order in anisotropy. The total phase and the total real exponent T of the 
influence functional are obtained by summing over k of ( |3.35| ) and ( |3.34| ) respectively. In 



performing these sums, divergent expressions will arise which will have to be regularized and 
renormalized. 

The above equations clearly show the history dependence of the influence functional since 
they are written in terms of time integrals of functions that depend on f3±{r]i) . Therefore, 
the phase and the real exponent are junctionals of the zero order WKB histories. 

Notice that since in this model we have more than one minisuperspace degree of freedom, 
even within a WKB branch, we have a whole family of trajectories rather than a single 
trajectory. So as far as the solutions of the Hamilton- Jacobi equation is concerned, this 
implies restricting ourselves to the family of trajectories given by the solution of ( 3.26) with 
a fixed value of b + /\b\. These are a family of parallel straight lines with the slope fixed by n 
and different (3 intercepts. We note that in the configuration space of the a — /3 + plane, one 
and only one trajectory passes through each point. Hence each point in configuration space 
can be associated with an entire history, and thus X is a functional of two histories as in the 
Brownian motion example. 

As it stands, T(r, r') is still not in a form that can be put in one-to-one correspondence 
with the !F(x, x') of the QBM problem, because the latter is an explicit function of time, 
whereas in the former, the WKB time is defined through ( |3.7| ) as a function of the coordinates 
r. The definition of 77= constant surfaces depends on the choice of the hypersurface in 
minisuperspace on which the initial condition of the wave function is specified. In our case 
the initial conditions were specified on a a = constant hypersurface. Thus our constant 
WKB time hypersurfaces are those with a= constant. Now, let us specialize to the situation 
where X(r, r') is evaluated on two points such that a = a'. From the above discussion then 
we know that this implies that rj = r)' '. The two histories, P±(rji) and P'±(r}i) that enter 
into the calculation of the influence functional are the parallel lines (with slope determined 
by (n)) , passing through the points (a,j3 + ) and (a,j3 + ') respectively. Now, the influence 
functional can be written as Z(/3±,P±, rj) and can finally be compared with that of the QBM 
problem. 
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4 Fluctuations in Quantum Fields and Dissipation of 
Spacetime Anisotropy 

4.1 Regularized Influence Action 



It has been pointed out in E3| that the influence action T is identical to the Schwinger- 
Keldysh (or Closed Time Path) effective action which is a functional of two histories and can 
be computed using diagrammatic techniques. Thus T is esentially the same as the quantity 



given by (3.11 ) in j4(J, with (3 and j3' corresponding to (3^ ,(3^ in the CTP context, where 
the + and — superscripts refer to the positive and negative contour branches respectively. 
This identification is useful as it connects with the well-known results in semiclassical gravity 
]40| . This connection provides both conceptual and technical advantages as it offers clearer 
physical interpretations of the results in quantum cosmology and makes available many 
results obtained previously in the application of the CTP formalism in quantum field theory 
in curved spacetimes. 



We now proceed to evaluate f and 9 by summing the equations (|3.34|) and (|3.35|) over 



all modes {k} subject to the restriction a = a'. Some of the mode sums appearing in these 
expressions are divergent and hence need to be regularized. The regularized influence action 



for this problem can be calculated using Feynman diagram |4(| or dimensional regularization 



techniques ||52|| . The phase of the influence functional can be written as 

e = v div + r ren (4.i) 

where Fdiv and T ren represent the divergent and finite contribution to the phase respectively. 
Fdiv (obtained as terms containing the 1/e factor in dimensional regularization , where e = 
n — 4 and n is the dimension of spacetime ) is given by |52| 



T dlv = J drndrniPij - P' tJ )(vihdUvi ~ V2W 3 + (4.2 

where 

r+00 

^/j.-..(rii — no) = I 

2tt 



-a; 4 



4(4nf(n 2 - 1) e 



(4.3) 



Tdiv can be rewritten as 



1 f •■ 2 •■ 2 

T div = — 727 o TT~ / d Vi[(3i ~P'i] + surface terms, (4.4) 

4(47r) [n z — lje J 

where the surface terms can be written as integrals of total derivatives of functions of (3 and 
(3' and can be discarded. As it stands this explicitly divergent term cannot be absorbed 
by renormalization of the bare coupling constants present in the original action since from 
(3.17) we see that no term of this higher derivative form appears there. Hence we follow 
the usual procedure used in quantum field theory in curved spacetime of first dimensionally 
regularizing the effective action, modifiying the original classical action by adding appropriate 
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counterterms to cancel the divergence, and finally taking the limit e — > 0. The modified 
classical action including the counterterms up to second order in j3 is given by |~7 



S 



drj 



-6Md 2 + 



1 f/^\ 4 /^\ 2 
{180(4tt) 2 } \ V ~ 3 V 



+ 



J dr] (m/?V + 



{180(4tt) 2 } 



Ze^p 2 + 3ln{iia)p - { {-)0 l + (-) 2 - (3 



a 2 o 2 



1.5) 



where \i has dimensions of mass and sets the renormalization scale. The total phase of the 
density matrix is now given by 



S{a,P)- S(a,/3') + Q(a,f3,P') 



(4.6) 



Inserting (4.5) for S in the above expression we notice that the pole term in e cancels exactly. 

. 

The rest of the exponent, T ren and T, is finite: 



and 



dvidrjzWii - (7/i)7(77i - ~ (3' iJ )(v2 



where the kernels 7 ren and 7 are given by 

1 



1ren(v) 



60(4tt) 



and 



1 



+ °° ^ e- ^ log(z^^ 
00 2n ji 



-°° du 



— 7TUJ COS LOT]. 
2,71 



(4.7) 
(4.8) 

(4.9) 
(4.10) 



60(4vr) 2 

Notice that the kernel 7(77) is even whereas jreniv) contains an odd and even part given by 



and 



'feveri 



1 



60(4vr) 2 Jo 
1 



60(4tt) 



-°° dw 
2tt 

00 duj 



iru smurj 



to 



— uj costur/ln — 

00 Z7T \X 



(4.11) 



(4.12) 



The kernel 7 ren (r/) is manifestly real and can also be seen to be causal [40 



Note that T and T ren play distinct roles here. T is responsible for the decoherence be- 
tween alternative histories /3 and j3' in the sense that it suppresses the contribution of widely 

1 However, we would like to add a cautionary note at this point. We are assuming without proof that the 
R 2 type terms can be added as counterterms at this level after making the WKB ansatz. Since addition of 
such terms at the level of the quantum cosmology Hamiltonian which was our starting point involves the 
introduction of new canonical degrees of freedom, the validity of this assumption is not entirely clear. 
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differing histories to the influence functional, and hence suppresses the off diagonal terms 
of the reduced density matrix . This feature and its connection to particle production was 
explored before in [52, 45]. On the other hand, when we attempt to derive the effective equa- 
tion of motion for (3 by varying the effective action S e ff, only T ren contributes to generating 
the equation of motion. The equation of motion obtained under such variation is identical 



to the real, causal dissipative equation for (3 obtained by Calzetta and Hu in [40|. In fact, as 



we will show more explicitly later, T ren provides the dissipative contribution to the equation 
of motion. Thus in the present form of the influence functional T ren contributes only to the 
equation of motion and not to decoherence, and T contributes only to decoherence, and not 
to the equation of motion. However, in the following we will show how T also plays the dual 
role of generating noise and will indeed contribute to the effective equations of motion with 
a stochastic source. 



4.2 Correspondence with QBM 

Now that we have the complete form of the influence functional, we can proceed to compare 
its exponent given by (|4.7|) and (|4.8|) with that of (|2.9| ) of the QBM problem. We can see 
that it corresponds to the k = 2, f(x) = x case in (2.9) with the identification /3j = x and 
q n = Xk- It is by no means obvious that our cosmological example should correspond to 
f(x) = x, i.e, the linear coupling case, because in our approximation we had retained up to 
quadratic terms in the anisotropy. In fact, from (|3.32|) we see that the system-environment 
coupling contains terms quadratic in f3 as well as a quadratic coupling in velocities, which is 
not even covered by our Brownian motion model. However, though these terms are originally 
present, when correctly dimensionally regularized, the terms proportional to A2 that contain 
the non-linear coupling vanish. Hence we are left with only an effective linear coupling in 
the anisotropy. The local potential renormalization terms AV's can be identified with 
in the cosmological case and we have already dealt with the renormalization. Using the time 
reflection symmetry of the kernel 7 we obtain 

n 771 

f = 2 JdrfiJ drfyjfa - m) - P nJ )(V2) (4.13) 


and for the phase T^ ren ) , using the time reflection properties of ^oddiv) an d leven(v) we can 
rewrite it as 

rn rv 

Tren = J J *7l*/2 Kj{Vl)l{Vl ~ fy)/? 3 (%) 

- J J ^1^2/3^(771)7(^-772)^(^2) 

+ 2 J d m J d V2 - pyXrn) lodd (Vi ~ %) + (3 ,IJ )(V2) (4.14) 


where 

7(771 - V2) = leven(Vl ~ ^) ~ 1 'add(Vl ~ V2) SQU^l ~ T] 2 ) (4.15) 
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is an even kernel. Now we must compare the expressions ( 4.13 ) for the real exponent and 



the phase (|4.14 ) with the corresponding expressions in the influence action ( [2.9|) for the 



Brownian motion case in order to properly identify the noise and dissipation contributions. 
Comparing the real exponents we see that the noise kernel for the anisotropy in this case is 
given by 

\ r+oo 

U{T1) = = %fi„ 27™ C ° S ^ (4 - 16) 

In trying to compare the phase terms we notice that the third term in ( |4.14| ) is indeed of 
the form of that in (2J3) and we can identify the dissipation kernel n(rj) for the cosmology 



case as 

//(Ty) = -2 lodd (r]) (4.17) 

and it is manifestly odd in time. 

The regularized influence action can therefore be written as 

r(/3, p) = f f d Vl dr, 2 faimmm - riPKm) 

V rv 

J dr} x dri 2 (7/1)7(7/1 - r] 2 )(3 m {r] 2 ) 
vi 

d Vl J d^fa-Tfe) (/^ + /?' l '% 2 ) 


V 



+ i / d Vl d V2 u(r ]l - V2 ) (^-/? w )(r/ 2 ) (4.18) 







The first two terms contribute a non-local potential to the effective action but do not con- 
tribute to the mixing of j3 and j3' histories like the third and fourth terms. We will now show 
in some greater detail that the third term with the kernel /1 that is odd in the time domain 
contributes to the dissipation and the last term containing v is associated with noise. 



4.3 Noise 

Let us first concentrate on the fourth term. Its contribution to the influence functional is 
given by 

exp[- f d Vl J d V2 u( Vl - r/ 2 ) (/F - I3 nj ){r l2 )} (4.19) 


We will proceed in exact analogy with the analysis of noise in the case of QBM described 
in Sec. 2.2. The term in (4.19) can be rewritten using functional Gaussian identity (2.11) 
which in this case states that the above expression is equal to 

v 

J D^V^expii J dieKi) (Pij ~ P'ijM)} (4.20) 
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where 

v v 

V[i] = P exp[- J drn J d^-CM^ivi ~ ^2)^(^2)] (4.21) 


is the functional distribution of £(77) and Po is a normalization factor given by 

v v 

Po' 1 = J Dt(v)exp[- J d Vl J d^M^ivi - V2)C j (V2)}. (4.22) 



The influence functional can then be written as 

e* r = jD^r])V[^expif[P,P'^} 

= <expif[(3,p',^>^ (4.23) 

where the angled brackets denote an average with respect to the stochastic distribution V[£\. 
The modified influence action r[/3, /?',£] is given by 

- f f d Vl drj 2 P'^mmvi - i»)/?*'(f») 



- fdr/Pfflfc • Jdr/FfflFy (4.24) 

The term coupling a stochastic source £ to (5 will manifest itself as the noise in the equation 
of motion derived from this effective action. We see that the influence action V can be 
written as an average of T over this stochastic distribution function. 

The reduced density matrix can thus also be written as a stochastic average 

Pred {f3,f3'}=<e l§ ^'^>t (4.25) 

where the full effective action S e ff is given by 

S eff = S[a,(3]+ jdr ] 'e ] (v')^-{S[a,P'}+ J drj'^Wij} 

V Vi 

J d Vl J dwfa - m) W ij + P' lJ )(V2) (4.26) 


and S is given by (4.5). Our relevant equations of motion will be derived by varying S e ff. 
From this equation we can view £(77) as an external stochastic force linearly coupled to f3, 
though the linearity is a feature specific to truncation of the perturbation series at quadratic 
order in the effective action. In general we will have non-linear coupling. 
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Since the distribution functional (|4.21|) is Gaussian, this is a Gaussian type noise, which 
as in (2.18) , is completely characterized by 

<£fa)> € = o 

<Z(m)£(m)> s = v(vi-V2) (4.27) 

Therefore the non-local kernel f(r]i — 772) is just the two-point time-correlation function of 
the external stochastic source Since this correlation function is non-local, this noise 

is colored. As suggested in ||2|, [H| we believe this is a rather general feature of noise of 
cosmological origin. 



4.4 Fluctuation-Dissipation Relation 

Now that we have identified the noise and dissipation kernels v{rj) and Joddiv) respectively, 
we can go ahead and write down the fluctuation-dissipation relation in analogy with the 



quantum Brownian model |26|, |57[] . Defining 



/i(r/) = -2lo d M = (4-28) 

The fluctuation-dissipation relation has the familiar form given by (2.30) 

00 

v ( v ) = [ dr}'K{r] - v'hiv') (4-29) 
where the FD kernel K{rj) is given by 







00 



K(v) = I — uj cos lot] (4.30) 

7T 







This supports the conjecture of [[G| that there exists a fluctuation-dissipation relation for 
the description of the backreaction effect of particle creation in cosmological spacetimes. 
We see that the FD kernel is identical with that given by fl2.36| ), which is given for more 



general system-bath couplings of the form (|2.3| ), but with the bath at T = 0. Hence this 
also vindicates the previous observation [^, ^ that the zero temperature fluctuation- 
dissipation relation is insensitive to the nature of the system-bath coupling. Since we have not 
taken the bath at a finite temperature, thermal fluctuations play no role in the above relation 
and it summarizes the effect solely of quantum fluctuations. Effect of thermal fluctuations 
can be included easily and we expect a FDR to hold for finite temperature particle creation 
and backreaction as well. 
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5 Particle Creation, Noise and Backreaction 



5.1 Particle Creation 

In this section we would like to examine in some detail the relationship between the noise 
and dissipation kernels and particle creation from the vacuum. We would also be interested 
in comparing this approach to that in and where the relationship between particle 
production and anisotropy dissipation was discussed in some depth. 

Let us first concentrate on the noise term . Since we know that the noise term comes 
from the real part V of the exponent of the influence functional, we will analyze this part 
and try to rewrite in a form such that it is easy to identify the part associated with particle 
production. It can be shown |52| that T(/3, (3') can be rewritten as 



m?) = I 7^st^\ b m - B k (P')\ 2 (5.1; 



d 3 k 1 
(2^4^' 
where 

B k {f3) = P d Vl — [2/%(r/i)A^] (5.2) 

and uj k = \k\ = (J2ih 2 ) 2 - As we may recall from (3.33), the term 2f3 i j(r]i)k l k : ' = Ai is the 
expansion of the natural frequency to the first anisotropy order. One can of course go to 
higher orders. 

Now we can show that a close relation exists between the B(f5) function and the Bogol- 
ubov coefficients associated with the particle creation that takes place as a consequence of 
the anisotropy evolution [17, 15]. This can be seen as follows. 

The conformally related massless scalar field X = a<3> in our model can be decomposed 
into modes as 

x= Iwf e '^ M (5 - 3) 

Xk satisfies the following equation to first order in anisotropy 



d 2 Xk 



+ (u k z + 2(3 l3 ( Vl )k l k J )xk = (5.4) 



drf 

The solution to the above equation ( again to first order in f3) is given by 



where 



-Xk in \v) fdrn —mji^kW} e 2 ™ (5.5) 



Xk in (v) = -nj=^ n (5-6) 



25 



is the solution to (|5.4 ) with = and corresponds to the 'in' conformal vacuum in the far 



past. Assuming that the anisotropy is switched off at time r\ the term on the left hand side 
of ( |5.5| ) can be associated with Xk out (v) the 'out' vacuum. As is well known, the "in" and 
"out" basis can be related in terms of Bogolubov coefficients a:*; and (3k as 

Xk° ut (v) = a kX k n (v) + PkXk m *(v) (5-7) 
Comparing (5.5 ) and (5.7 ) we can identify the (3k Bogolubov coefficient as 

k = fd Vl J-^fa)^') (5.8) 

As we see from its definition, the function B is proportional to this Bogolubov coefficient 

$• 

Bk((3) = u k e- 2 ^$k (5.9) 

Thus 

m?) = J 77^34^ -^i 2 » ( 5 - 10 ) 

where (3k and (3' k are the Bogolubov coefficients associated with the anisotropy histories (3^ 
and respectively. It is obvious from (5.19 ) that the noise will be non-zero only provided 
(3k 7^ (3'k, i-e, if there is different amounts of particle production along the two histories. 
Since this term is also associated with decoherence this is also a necessary condition for 
decoherence to occur. This has also been noticed from a slightly different point of view in 



IH |9 



This demonstrates a connection between the process of particle production and the noise 
or fluctuation. 



5.2 Einstein-Langevin Equation 



We will now show how this noise can be incorporated into the equation of motion as a 
Langevin type equation. In this process we will also demonstrate the role of the kernel \x in 



providing dissipation. The key difference from the earlier treatment J4C§ is that the equation 
of motion will be derived from the quantity S e ff((3, (3', £) rather than the "noise averaged" 
quantity S e ff((3, (3'). This has also been discussed in other contexts in |26|, [45], |6[]. The first 
step is to write S e ff((3, /?',£) in terms of the following variables 



Pij — 

A = 

The equation of motion is then derived as 



fly - &i 



13) 



A) 

SA 



A=0 



(5.11) 



(5.12) 



26 



yielding 



d 



2M—(a 2 (3 ij ) + 



1 



d 2 



]Pijln(fla)} + 



d 



<V " lJ> ' 30(4vr) 2 ^ 2L 3 ^ n 90(4vr) 2 rfr/ 
/ dr}xlren{r] - m)Pij(rn) = -jij(r]) + £ij(rj) 



/CL\ /CL 
(a) + (n 



(5.13) 



Here j%j is an external source term added in order to switch on the anisotropy in the distant 

where similar equations were 
we find that they 



past fl7| . It is worth comparing these results with those in 
deduced from the CTP effective action. Comparing ( J5.13|) with (3.18) in 
are exactly the same except for the stochastic force £y on the right hand side. The real and 
causal kernel K4 there (including the numerical factor l/[30(47r)] 2 ) is identical to our kernel 
7 ren . We will show that the odd part of this kernel can be associated with dissipation. One 
could in fact interpret (3.18) obtained by Calzetta and Hu as (|5.13| ) averaged with respect 
to the noise distribution. Since this is a Gaussian noise, < £ >= 0, we obtain (3.18) of [i0|, 
where the /3's are also to be interpreted as noise-averaged variables. In this sense, we have 
gone beyond previous analysis in extracting the underlying stochastic behavior that is lost 
in the smoothed out average version given in |40] . 



To make the analogy with a Langevin equation more explicit it is convenient to integrate 
( |5.13j ) once with respect to 77. This gives the following equation 



-2Ma 2 f3ij + 



d [(3ijln(/j,a)] + 



1 



+ 



30(4tt) 2 drj 
dr) 2 / drji 



90(4tt) 



(-) +(-) 



(5.14) 



where ^(77) = -Jdrfj^rf) and ._Sij{r}) = fdrf^j(rf). 

Defining the variable = dfiij/drj we can write the above equation in the following form 



(M J 
drj drj 



where 



M 
k 

JCqij 



1 



30(4vr) 
-2Ma 2 + 



drj 



In(jla) 



90(4vr 
drj 2 I difrfto - r)i) 



fd\ ('a 

(a) + L 



dqjj 
drj! 



(5.15) 

(5.16) 
(5.17) 
(5.18) 



and d 2 f(rj)/drj 2 = j ren . This equation is identical in form to the equation (3.15) in J42| 
except for the term on the right hand side, which is indeed the stochastic contribution 
from the noise anticipated there. This equation is a generalized Einstein equation in the 
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Langevin form, in that there is a dissipative term in the dynamics and a noise term in the 
source. It has been conjectured |42| and shown [(45| that in a more complete description of 



semiclassical gravity the semiclassical Einstein equation driven by the expectation values of 
the energy-momentum tensor should be replaced by an Einstein-Langevin equation, where 
there is an additional stochastic source arising from the fluctuations of quantum fields. The 
conventional semiclassical Einstein equation is in this sense, the simplified mean-field theory. 

5.3 Dissipation and Backreaction 

Equation ( [5.15|) is in the form of a generalized damped harmonic oscillator driven by a 
stochastic force Sij. (Of course the generalized mass M and spring constant k are time 
dependent, so strictly speaking it has the damped harmonic oscillator analogy only when 
these quantities are positive, as was also pointed out in |6(J.) 



The second term on the left hand side of ( |5.15 ) represents the damping term involving a 



non-local (velocity dependent) friction force. That this term is associated with dissipation 
can be quickly seen as follows . In the Fourier transformed version of a damped harmonic 
oscillator equation the imaginary term is associated with dissipation. Writing (|5.15|) in 
terms of the Fourier transform qij{uj) = J drje~' lu '' n qij(r]) we notice that the only imaginary 
contribution comes from the second term on the left hand side, which can be written as 

J ATI UJ £ 

where 7 re n(^) is the Fourier transform of Jreniv) defined in ( fOQ . Thus we see that the 
dissipation is associated with the imaginary part of ^ ren {uj) or equivalently with the odd 
part of the kernel ■jreniv) given by j dd defined in ([4.11]) . This is consistent with our earlier 
identification of ^ dd as the dissipation kernel from the form of the influence functional 



compared with the Brownian motion case. In fact, as in [4C , 41] we can isolate the generalized 
(frequency dependent) viscosity function (,{oj) by writing 

i((Lo)Loqij(uj) = i\m.^ Ten (uj)q ij (uj) (5.20) 

From ( |4.9| ) we can identify (,{&) as 

« w > - <5 - 2i) 



which is identical to that found in [42] and which is not surprising, since our kernel 7( ren ) in 
( |5.13| ) and in (3.18) in |I(J are identical up to numerical factors. 

Once having made the identification of the velocity dependent viscous force in the equa- 
tion of motion , we can calculate the dissipated energy density by integrating F.v (with 
v = <jij acting as velocity) over all frequencies and come up with an expression identical to 
(3.18) in g| H. 

Pdissipation = [ ^[^AjM1[C(^VAjM]- (5-22) 
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which has been shown there to be identical to the total energy of particle pairs created 
by a given anisotropy history f3. In this way, we can see the connection between particle 
production and the dissipation kernel and hence the process of dissipation itself. Earlier 
in this section we had demonstrated the connection between particle production and the 
noise or fluctuation term. On the other hand, ( |4.29| ), the fluctuation-dissipation relation, 
embodies a relationship between the processes of fluctuation and dissipation of anisotropy. 
So this completes the full circle of connections among these processes. In a way, one can 
say that as a physical process, particle production is contributing to both the noise and 
dissipation, and of course these are two different manifestations of the loss of information 
due to integrating over the field modes. 



6 Discussion 

In closing, we would like to discuss the meaning of the FDR in semiclassical cosmology in a 
broader context and mention some related problems for future investigation. 

1) FDR under Finite Temperature and N on- Equilibrium Conditions 
In this paper we have discussed in detail the FDR in semiclassical cosmology under a 
zero temperature bath. A similar relation between the noise and dissipation kernels exists 
for baths at finite temperature. The form will be similar to that derived for the QBM 
problem in Sec. 2 p6 |. One can take the finite temperature calculation via the CTP for- 



malism and perform a similar analysis as we have done for the vacuum case and obtain 
the results explicitly. In reality both vacuum and thermal bath results will enter into the 
picture,^ since once particle creation commences, given sufficient time and assuming some 
(collisional) interaction amongst the created particles, the bath will soon acquire a finite 
temperature character .0 This heat-up process is expected to happen quickly near the Planck 
time, especially so for anisotropic universes, as particles are created profusely there |12], |10j, 
generating a large amount of entropy p|l ^9], ^2], [F(], [nj . The copious creation of particles 



near the Planck time is accompanied by large fluctuations and noise, and it induces a strong 



backreaction on the spacetime dynamics, dissipating the anisotropy rapidly |L2|, [15], [L7[]. The 
weaker anisotropy in the universe's expansion induces lesser particle creation. The lower 
particle creation rate is accompanied by a smaller fluctuation and noise, which in turn gives 
weaker dissipation of spacetime anisotropy. The surviving anisotropy would continue to sus- 
tain particle creation, albeit in much smaller amounts. And this goes on. (The backreaction 

2 As has been discussed earlier plj , the energy density of the quantum field at any moment will contain two 
parts. There is a zero temperature component and a finite temperature component, the former corresponds 
to spontaneous creation from the vacuum, and the latter is of the nature of stimulated creation from particles 
already present [ffl |62| p3|, |64j. 

3 A finite temperature bath at every moment is only an idealization. To use a finite temperature description 
one has to discriminate the conditions for the bath to thermalizc, and for the system to be equilibrated with 
it. These vary with the nature of the bath (massive or massless, linear or nonlinear interactions, spectral 
density) and the form of interaction between the system and the bath. See the analysis of |35| ^36] ^37j 
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follows a Lenz law behavior which was expounded in earlier studies ]7J, [72], 59 1.) At each 
stage we expect to see a balance between the rate of particle creation and the strength of 
fluctuations and dissipation. 

In reality the spacetime-field combined system involving particle creation exists in a 
highly non-equilibrium state. To give a quantitative description of the above processes 
one needs to describe the dynamics of the actual statistical state of both the system and the 
environment in a self-consistent manner, which is a highly non-trivial problem. What we have 
described in this paper is only the first step, which depicts the effect of particle production 
from a vacuum (zero-temperature) bath. The interaction of created particles and how they 
alter the environment (e.g., thermalization) is not accounted for. In the second step one needs 
to also examine the evolution of the environment (quantum field) taking into consideration 
the effects of spontaneous and induced particle creation, their interaction and the entropy 
generation processes, all in the context of a changing background spacetime whose dynamics 
at each moment affects and is also affected by the activities of this environment. 

Despite all these complexities, even in highly nonequilibrium conditions we expect that 
a generalized FDR (in the form given in this paper) will still hold and be useful to guide us 
on understanding the complex physical processes in the system and the environment. From 
the above depiction of the physical scenario and from previous studies of the statistical 
mechanics of quantum field processes in cosmology, one can see that there exists a balance 
between particle creation (in the field) and its backreaction (on spacetime), which can be 
attributed to the interlocked relation between fluctuations and dissipation. There is also 
a mathematical justification: it is a relation between the real and imaginary parts of the 
effective action for the open system. Similar in nature to the optical theorem in scattering 
theory or the Kramers-Kroning relation in many body theory [fi5| |, these relations describe 
the dissipative and reactive parts of the response function of an open system to influences 
from the environment. They are of a categorical nature because they originate ultimately 
from the unitarity condition of the dynamics of the combined closed system. They only take 
the form of dissipation in the open system because we have identified a certain subsystem 
as the system of interest and decided to follow its effective dyanamics; and they take the 
form of fluctuations in the environment because we refer to them in reference to the mean 
value of the environment variables, the remaining information is downgraded in the form 
of fluctuations. Had we decided not to coarse-grain the environment, or choose to observe 
the two subsystems with equal interest and accuracy, such a relation governing the mutual 
reaction would still exist, except that the concepts of dissipation and fluctuations will no 
longer be appropriate. (Both subsystems will be governed by equations of motion in the 
form of an integral differential equation, and treated in a nondiscriminate and balanced way. 
See, e.g., [[fj). 

In the context of semiclassical cosmology, the open system is the spacetime sector, whose 
dynamics is influenced by the matter fields. The expansion of the universe amplifies the 
vacuum fluctuations of the matter field into particles, which act as the source in the Einstein 
equation driving the universe. The averaged effect of particles created imparts a dissipative 
component in the spacetime dynamics, and the fluctuations in particle creation constitute 
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the noise. The particular forms of the dissipation and noise kernels and their effects may 
vary under different conditions- zero or finite temperature, equilibrium or non-equilibrium- 
but the existence of such a relation between the fluctuations in the matter field and the dis- 
sipative effect on the spacetime dynamics should remain. We will have opportunities later 
to explore related problems which can shed more light on these issues. 



2) Relation with FDR in Spacetimes with Event Horizons 

As we mentioned in the Introduction, our search for a FDR in cosmological spacetimes 
without event horizons was inspired by Sciama's proposal to view the Hawking and Unruh 
effects as manifestations of a fluctuation- dissipation relation between the field quanta and 
the detector response. De Sitter universe is an important class of cosmological spacetimes 
with event horizons. For this one can use the thermal property of the field to perform a 



linear response theory (LRT) analysis for the derivation of the FDR [37|. Our derivation 



here based on the influence functional formalism attacks the problem at a more basic level, 
where equilibrium condition between the system and the environment is not necessarily 
present at every stage. It is of interest to compare the results between the equilibrium 
limit of the IF or the CTP formalisms and that of LRT. This can be done explicitly by 
carrying out an analysis similar to this paper on the de Sitter universe and see how the 
FDR obtained from the IF compare with that from the LRT. Formally this would render 
explicit the relation between the IF formalism to (non-equilibrium) statistical field theory 
and perturbative thermal (finite-temperature) field theory. 



More meaningfully, as was originally concieved by one of us [|74], f|2|, this would pro 



vide a channel to generalize the conventional way of treating Hawking effect associated with 
black holes and accelerated observers based on thermal propagators and event horizons to 
non- stationary conditions. This involves cases like non-uniformly accelerating observers and 
realistic collapse dynamics, where an event horizon or Euclidean section does not always 
exist but is dynamically generated. Our motivation for finding a way to treat these more 
general conditions is to seek a deeper meaning to the Hawking effect, and through it to 
explore the subtle connection between quantum field theory, relativity theory and statistical 
mechanics. In our view, the open system concept explicated by the influence functional 
formalism provides a more solid basis to understand its statistical mechanical meanings and 
a broader framework to tackle the less unique situations which cannot easily be treated by 
purely geometric means, powerful and elegant as they are. It also brings the effects of quan- 
tum fields on observer kinematics and spacetime dynamics in line with the more common 
statistical mechanical phenomena involving ordinary matter. These problems are currently 
under investigation. 



3) Related Problems in Semiclassical Cosmology and Inflationary Universe 
For particle creation-backreaction problems similar to the Bianchi-I model studied here, 
one can obtain similar results for other matter fields in other types of spacetimes of astro- 
physical or cosmological interest. An example is a massless minimally-coupled scalar field in 
a Robertson- Walker or de Sitter universe. It mimics the linearized graviton modes and has 
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practical use for the description of primordial stochastic gravitons. The particle production 
problem was first studied by Grishchuk |T]| , the backreaction by Grishchuk and Hu and 
Parker [14 1 via canonical quantization methods, and by Hartle [18] and Calzetta and Hu 
|40fl via the in-out and in-in effective action method respectively. The influence functional 
approach expounded here would enable one to get from first principle the entropy generation 



from graviton production |70|, 71], and the noise associated with them, which is related to 
the fluctuations in graviton number |64 , |45| , On the aspect of backreaction in semiclassi- 
cal cosmology, one can also derive the Einstein-Langevin equation for the study of graviton 
production and metric fluctuations. This problem is currently pursued by Calzetta and Hu 

m 

A related problem of interest is the evolution of the homogeneous mode of the inflaton 
which describes the inflation mechanism |76| and the inhomogeneous modes as progenator 
of structures in the early universe. The influence functional method was used by Hu, Paz 
and Zhang ]57J] Laflamme and Matacz [|77[] and others to discuss the decoherence of the long- 
wavelength sectors of the inflaton, and the origin of quantum fluctuations as noise for the 
galaxy formation problem. Our result here provides an example for the consistent treatment 
of the evolution of these modes, their intereaction, and their backreaction on the spacetime, 
which can offer some physical insight into the no-hair type of theorems in inflationary uni- 
verse. These problems are under study by Matacz, Raval and the authors. 

4) Minisuperspace in Quantum Cosmology as an Open System: Geometrodynamic Noise 
and Gravitational Entropy 

We have discussed the question of the validity of the minisuperspace approximation [|5B] 
in quantum cosmology 179], |80|] , wherein only the homogeneous cosmologies are quantized 



and the inhomogeneous cosmologies ignored [^T|. We used an interacting quantum field 
model and calculated the effect of the inhomogeneous modes on the homogeneous mode via 
the CTP effective action. This effect manifests in the effective equation of motion for the sys- 
tem as a dissipative term. For quantum cosmology, this backreaction turns the Wheeler-De 
Witt |82] equation for the full superspace into an effective equation for the minisuperspace 
with dissipation. Extending the CTP to the IF formalism as is done here, one can derive the 
noise associated with the truncated inhomogeneous cosmological modes. One can also define 
an entropy function from the reduced density matrices, which measures the information loss 
in the minisuperspace truncation. These can perhaps be called geometrodynamic noise and 
gravitational entropy. It would be interesting to compare this statistical mechanical defini- 



tion with the definition suggested by Penrose [S3 in classical general relativity and by one 
of us in the semiclassical context f72|| . Some initial thoughts on this problem are described 
~~4j, while details are to be found in 



in 
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